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在 解释 物种 丰富 度 
昌 蛇 类 物种 丰富 度 的 地 理 分 布 格 
该 研究 采用 主 节 导 基于 特征 值 的 空域 数据 ， 并 同 生态 
丰富 度 在 经 、 纬 度 上 呈现 多 峰 分 布 格局 : 高 的 
于 青藏 高 原 、 北方 E Jn 353 、 TEF Jr HHY 
FE 富 度 变 化 ， 分 忆 蛇 类 物种 丰富 度 格局 的 3 向 因子 是 归 
。(3) 模型 选择 结果 显示 元 回归 分 析 中 ，P<0.05 fi -化 植被 
的 模型 是 解释 蛇 类 物种 的 最 优 模型 。 这 说 度 
甘于 中 国 蛇 类 物种 丰富 成 因 研究 的 复杂 性 , - 
择 和 人 类 活动 的 景 间 尺 度 上 对 蛇 类 物种 
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进 果 表 明 : (1) 中 
位 于 东洋 界 亚热带 、 热 带 ， 丰富 度 较 1 
等 ; (2) 多 元 回归 分 析 能 解释 56.5% 的 蛇 
化 植被 指数 、 最 冷 季 降 水 量 和 年 温差 
指数 、 最 冷 季 降水 量 和 年 温差 ) 组 成 
不 同 生 态 因 子 共同 作用 的 结果 。 
P， 需 重视 各 假说 中 影响 因子 的 选 
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Geographic patterns and ecological factors correlates of 
snake species richness in China 


CAI Bo?, HUANG Yong ?, CHEN Yue-Ying !, HU Jun-Hua !, GUO Xian-Guang !, 
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(1. Chengdu Institute of Biology, the Chinese Academy of Science, Chengdu 


100049, China; 3. Guangxi Botanical Garden of Medicinal Plants, Nanning 530023, China) 


610041, China; 2. Graduate School of Chinese Academy of Science, Beijing 


Abstract: Understanding large-scale geographic patterns of species richness as well its underlying mechanisms are 
among the most significant objectives of macroecology and biogeography. The ecological hypothesis is one of the most 
accepted explanations of this mechanism. Here, we studied the geographic patterns of snakes and investigated the 
relationships between species richness and ecological factors in China at a spatial resolution of 100 kmx100 km. We 
obtained the eigenvector-based spatial filters by Principal Coordinates Neighbor Matrices, and then analyzed ecological 
factors by multiple regression analysis. The results indicated several things: (1) species richness of snakes showed 
multi-peak patterns along both the latitudinal and longitudinal gradient. The areas of highest richness of snake are tropics 
and subtropical areas of Oriental realm in China while the areas of lowest richness are Qinghai-Tibet Plateau, the 
grasslands and deserts in northern China, Yangtze-Huai Plain, Two-lake Plain, and the Poyang-lake Plain; (2) results of 
multiple regression analysis explained a total of 56.596 variance in snake richness. Among ecological factors used to 
explore the species richness patterns, we found the best factors were the normalized difference vegetation index, 
precipitation in the coldest quarter and temperature annual range ; (3) our results indicated that the model based on the 
significant variables that (P«0.05) uses a combination of normalized difference vegetation index, precipitation of coldest 
quarterand temperature annual range is the most parsimonious model for explaining the mechanism of snake richness in 
China. This finding demonstrates that different ecological factors work together to affect the geographic distribution of 
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snakes in China. Studying the mechanisms that underlie these geographic patterns are complex, so we must carefully 
consider the choice of impact-factors and the influence of human activities. 


Key words: Macroecology; Ecological hypothesis; Biogeography; Spatial distribution patterns; Energy hypothesis; 


Habitat heterogeneity hypothesis 




















物种 丰富 度 大 尺度 地 理 格局 及 其 成 因 是 生态 
学 和 生物 地 理学 研究 的 重要 目标 之 一 (Gaston， 
2000)。 目 前 已 有 多 种 假说 解释 物种 丰富 度 地 理 格 局 
的 形成 机 制 (Currie, 1991; Gaston, 2000; Hawkins et 
al, 2003a,b; Rahbek, 2005)， 其 中 ， 受 到 广泛 关注 的 
生态 假说 (ecological hypothesis) 主 要 包括 生产 力 假 
说 (productivity hypothesis)、 环 境 能 量 假说 (ambient 















































energy hypothesis). 3E7? lij f Ut (freezing tolerance 
hypothesis)、 水 能 动态 假说 (water-energy dynamic 
hypothesis) 及 生境 异 质 性 假说 (habitat heterogeneity 
hypothesis)(Gaston, 2000; Terribile et al, 2009)。 生 产 
力 假说 认为 能 量 的 增加 提高 了 一 个 地 区 的 净 初 级 
生产 力 ， 导致 流 向 食物 网 各 营养 级 的 能 量 增加 , D 
证 了 更 多 个 体 的 存活 ， 从 而 提高 了 物种 多 样 性 
(Currie et al, 2004; Hawkins et al, 2003a, b; Wang et 
al, 2009; Wright, 1983); 环境 能 量 假说 认为 环境 能 
量 直 接 影响 生物 个 体 的 生理 调节 机 制 ， 从 而 影响 种 
群 增长 进而 影响 物种 丰富 度 ( Currie, 1991; Huang et 
al, 2011; Turner, 2004; Turner et al, 1987; Wang et al, 
2009); 寒冷 忍耐 假说 认为 多 数 物种 因 受 最 冷 季 低 
温 忍 受 力 影响 而 不 能 在 更 冷 的 地 区 生存 (Currie et al, 
2004; Hawkins, 2001; Hawkins et al, 2003a; Wang et 
al, 2009); 水 能 动态 假说 基于 水 分 对 生物 能 量 利用 
过 程 的 影响 ,认为 生物 多 样 性 的 大 尺度 分 布 格局 | 
水 能 动态 决定 ( Hawkins et al, 2003a; O'Brien, 1993; 
Wang et al, 2009); 生境 异 质 性 假说 认为 物理 因子 或 
生物 因子 的 复杂 性 孕育 较 多 生态 位 , 环境 越 复杂 ， 
空间 异 质 性 程度 越 高 ,， 小 生境 类 型 越 多样 ， 所 能 
持 的 物种 数 也 越 多 ( Cox & Moore, 2010 ; Jiang et al, 
1997; Kerr & Packer, 1997; Pianka, 1966)。 但 越 来 越 
多 的 研究 表明 物种 丰富 度 地 理 格局 的 形成 机 制 征 
复杂 , 不 能 仅 由 某 一 个 假说 来 解释 (Currie, 1991; 
Hawkins et al, 2003a; Rahbek, 2005)。 




































































































































































1999), 截止 2012 年 2 月 ,我 国境 内 记录 蛇 类 一 222 
种 (Gong et al, 2010; Guo et al, 2008; Guo & Wang, 
2011; Guo et al, 2011; Huang et al, 2012; Jiang & 
Zhao, 2009; Luo et al, 2010; Malhotra et al, 2011; Rao 
& Zhao, 2005; Vogel & David, 2010; Vogel & Luo, 
2011; Yang et al, 2011a, b; Zhang et al, 2011a; Zhang 
et al, 2011b; Zhao, 2006; Zhao & Adler, 1993), ZJ i 
全 球 蛇 类 物种 数 的 796 (http://www.reptile 
-database.org)。 关 于 中 国 蛇 类 物种 多 样 性 , 早 在 三 
年 前 的 《山海 经 》 中 就 有 记载 ， 其 分 布 格局 与 生态 
因子 关系 则 在 20 世纪 80 年 代 开 始 受 到 关注 (Shen， 
1985)。 尽管 在 一 些 两 栖息 行动 物 的 分 布 格局 研究 中 
涉及 蛇 类 , 但 尚 无 针对 整个 中 国 区 域 蛇 类 分 布 格局 
的 相关 研究 ,， 这 与 该 类 群 所 具有 的 独特 进化 地 位 很 
不 相称 。 

不 同类 群 物种 具有 不 同 的 生理 特性 , 适应 不 同 
的 生态 环境 ,其 分 布 格局 及 其 成 因 也 有 所 不 同 
(Gaston, 2000; Lin et al, 2009), 特别 在 全 球 气 候 变 
化 背景 下 ,更 需要 了 解 某 一 类 群 的 物种 多 样 性 分 布 
格局 与 生态 因子 关系 机 制 ， 以 利于 生物 多 样 性 保 
护 。 因 此 ， 本 文 将 探讨 中 国 蛇 类 物种 丰富 度 分 布 格 
局 及 其 影响 因素 , 结合 生态 假说 分 析 其 成 因 , 为 蛇 
类 多 样 性 保护 提供 科学 依据 。 


1 材料 和 方法 


1.1 物种 数据 

EA FÉ fr zy 99 X d 库 (http//www.reptile- 
database.org) 和 《中 国 蛇 类 》(Zhao, 2006) 的 分 类 系统 
为 依据 ， 本 研究 收集 中 国 蛇 类 物种 信息 , 包括 采集 
地 点 经 纬度 、 生 境 与 海拔 等 。 数 据 主要 来 源 于 中 国 
TEE LPS RETI) As [Museum of Herpetology, 
Chengdu Institute of Biology (CIB), Chinese Academy 
of Sciences], LAE CP Ed JG 8 66 H ENTE. H «€ 














































































































































































































































































































中 国 疆 域 广阔 ， 跨 越 古 北 与 东洋 两 大 动物 区 




















国 蛇 类 》 和 在 以 下 杂志 公开 发 表 的 文献 ， 如 











系 。 第 三 纪 晚 期 以 来 发 生 的 地 质 事件 , 特别 是 印度 
板块 北 移 与 青藏 高 原 隆起 ， 使 得 区 域内 地 理 环 境 复 
杂 , 气候 多 样 ， 物 种 (包括 蛇 类 ) 多 样 性 丰富 (Zhang， 









































(Zootaxa) . 《Bonn Zoological Bulletin) ~ 《Asian 
Herpetological Research》、《 两 栖息 行动 物 学 报 》、 
《动物 学 研究 》、 《四 川 动物 》、《 动 物 分 类 学 报 》 等 。 
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若 资 料 中 记录 有 经 、 纬 度数 据 ， 则 采用 原 资 料 的 经 、 
纬度 ; 若 无 ， 则 根据 资料 记载 地 点 (至 少 到 县 级 行政 
区 划 ) 使 用 Google Earth 查找 物种 分 布 的 经 、 纬 度 。 
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另外 , 我 们 还 邀请 相关 专家 对 原始 数据 审阅 、 更 新 和 
修正 。 数 据 整 理 截 止 2011 年 8 月 ， 共 获得 19 575 条 
记录 (图 D, 最 终 构 建 了 中 国 219 种 蛇 类 分 布 数据 库 。 
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图 1 基于 标本 馆 记 录 和 各 种 文献 收集 的 中 国 


Raw sample point of snakes based on museum collections and published literature 


Fig. 1 

















蛇 类 分 布点 图 





地 图 投影 基于 Alberts Equal-Area Conic Projection ( 双 标 准 纬 线 : N25° 和 N47°,， 中央 经 度 E105°)。 
The map used an Alberts Equal-Area Conic Projection (at the standard parallels of N25°and N47°, central meridian of E105°). 


12 生态 数据 

本 研究 共 采 用 31 个 生态 因子 作为 自 变量 ( 表 1), 
并 基于 生态 假说 将 其 分 为 5 组: 生产 力 假说 、 环 境 
Eb 量 假说 、 寒 冷 忍 耐 假说、 水 能 动态 假说 和 生境 异 
质 性 假说 (Currie, 1991; Hawkins et al, 2003a; Kerr & 
Packer, 1997; Terribile et al, 2009; Wang et al, 2009). 
相关 假说 和 生态 因子 如 下 : 

生产 力 假说 一 般 用 实际 燕 散 量 (mean annual 
actual evapotranspiration, AET) (Currie & Paquin, 
1987; Huang et al, 2011) 和 归 一 化 植被 指数 
index, NDVI) 
(Evans et al, 2006; Huang et al, 2011)fffj & . AET 的 空 
间 分 辨 率 为 30， 数 据 来 源 于 http:/www.grid.unep.ch/ 
data/grid/gnv183.html. NDVI 数据 来 自 1998—2008 
E, 空间 分 辩 率 为 1km, 来源 于 http;//westdc. 
westgis.ac.cno 


环境 能 量 假说 采用 温度 及 相关 指标 衡量 (Currie， 
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(normalized difference vegetation 





























1991; Huang et al, 2011; Turner et al, 1987; Wang et 
al, 2009). Æ 3C 3 H] E 235] il (°C (annual mean 
temperature, Tnem)、 等 温 性 (C) (isothermality, ISO). 
温度 季节 性 变化 (C) (temperature seasonality, TS) 、 

风速 (m/s)(mean annual wind, WIND) fj tg, t P 
Tue» ISO. TS 数据 来 自 1961—1990 年 ,空间 分 辩 
率 为 30"( 一 1 km2， 来 源 于 http:/www.worldclim. 
org. WIND 来 自 1981—1996 年 ， 空 间 分 辨 率 为 500 
m2， 来 源 于 http;//www.geodata.cn/Portal 。 
寒冷 忍耐 假说 常用 最 冷 季 均 温 、 最 冷 月 均 温 或 
年 最 低温 等 极端 低温 相关 指标 来 衡量 (Hawkins， 
2001; Hawkins et al, 2003a; Wang et al, 2009)， 又 因 
物种 对 低温 的 耐 受 能 力 受 所 处 环境 温度 及 其 变化 
影响 ， 因 此 ,同时 也 需要 采用 一 般 环 境 温 度 、 温 差 
指标 来 衡量 (Ricklefs & Schluter, 1993; Wiens & 
Donoghue, 2004)。 本 文采 用 年 最 高 温 (C)(mean 
annual highest temperature, Tmax)、 年 最 低温 ('C) 
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R1 表征 各 个 假说 的 生态 因子 


Tab.1 Tested ecological factors associated with the five hypotheses 





生态 假说 Ecological hypothesis 





生态 因子 Ecological factors 








生产 力 假说 

Productivity hypothesis 
环境 能 量 假说 5 
Ambient energy hypothesis 





寒冷 忍耐 假说 
Freezing tolerance hypothesis 

















Habitat heterogeneity hypothesis 


(mean annual lowest temperature, Tnin)、 最 冷 月 均 温 
('C)(temperature of coldest month, TCMmean)、 最 热 月 
均 温 (CC)(temperature of warmest month, TWM mean) ~ 
fcil] ECC (mean temperature of wettest quarter, 
TWetQmean) 、 最 旱季 均 温 ('C )(mean temperature of 
driest quarter, TDQmean) 、 最 冷 季 均 温 ('C )(mean 
temperature of coldest quarter, TCQmean)、 最 热 季 均 温 
( C )(mean 
7TWarmQaeam)、 人 年 温差 (CC )(temperature annual range, 
TAR), Hi 2:('C)[mean diurnal range-mean of 
monthly (max temp-min temp), TDR] 衡 量 ， 这 些 数 据 
来 自 1961—1990 4E, TEDIN 30" (— 1 km?) 
来 源 于 http://www.worldclim.org » 

水 能 动态 假说 是 基于 水 分 对 生物 能 量 利 用 过 
程 的 影响 ,认为 生物 多 样 性 大 尺度 分 布 格局 由 水 能 
动态 决定 , 通常 用 蒸 散 量 、 温 度 、 降 水 量 和 湿度 相 
关 指 标 来 衡量 (Hawkins et al, 2003a; Wang et al, 
2009). Zk SC HI 4E (C )(mean annual temperature, 
Tmean) SE 2E HUE (mm/y)(actual evapotranspiration, 
AET), WEZE (mm/a)(potential evapotranspiration, 
PET). fF [E 7K Œ (mm/a)(mean annual precipitation, 
AP). $ H D RS Œ (mm/a)(precipitation of wettest 
month, PWM), i^ H FER = (mm/a) (precipitation of 
driest month, PDM) ~ Æ Ti TE K 7K Œ (mm/a) 
[precipitation seasonality (coefficient of variation), 
PS]. iH Æ K k Œ (mm/a) (precipitation of wettest 
quarter, PWetQ)、 最 旱季 降水 量 (mmya) (precipitation 
of driest quarter, PDQ) . E 7$ z& KK 
(mm/a)(precipitation of coldest quarter, PCQ)、 最 热 季 




































































temperature of warmest quarter, 























































































































实际 蒸 散 量 (mm/y)E7、 归 一 化 植被 指数 (VDVD 


年 均 温 ('C)(Tmean)、 等 温 性 ISO)、 温 度 季 节 性 变化 (TS)、 风 速 (m/s)(WIND) 


年 最 高 温 ('C)(Tmax)、 年 最 低温 ('C)(Tmin)、 最 冷 月 均 温 ('C) (TCMmean)、 最 热 月 均 温 ('C) (TWMmean)、 最 潮 
季 均 温 ('C) (TWetQmeam)、 最 旱季 均 温 ('C) (TDQmewm)、 最 冷 季 均 温 ('C) (TCQua). BA ES 
CC) TWarmQmean)、 年 温差 ('C) (TAR)、 


年 均 温 (C)(Tiwean)、 实 际 蒸 散 量 (mm/a)(A4E7T)、 潜 在 共 散 量 (mm/a)(PET)、 年 降水 量 (mm/a)(4P)、 最 潮 月 降 














温差 CC) (TDR) 








水 能 动态 假说 雨量 (mmy/a)(PWM)、 最 旱 月 降 南 量 (mm/a) (PDM)、 季 节 性 降水 量 (mm/a) (PS) 、 最 潮 季 降水 量 

Water-energy dynamic hypothesis (mm/a)(PWetQ)、 最 早 季 降水 量 (mm/a)(PDQ)、 最 冷 季 降水 量 (mm/a)(PCQ)、 最 热 季 降水 量 (mm/a) 
(PWarmQ)、 相 对 湿度 (%)(REH)、 潮 湿 天 数 (mm/a)(WD) 

生境 异 质 性 假说 


土地 有 覆盖 类 型 数 (LUCC)、 平 均 海拔 (ELEmean)、 最 高 海拔 (ELEmsx)、 最 低 海拔 (ELEnmin) 


降 水 量 (mm/a)(precipitation of warmest quarter, 
PWarmQ)、 相 对 温度 (%)(relative humidity, REH). 
淹 湿 天 数 (mm/a)(wet-days, Number days with>0.1 
mm rain per month, WD) 衡 量 。 其 中 PET 的 空间 分 辨 
X 为 30, X W 于 http:/www.grid.unep.ch/data/ 
grid/gnv183.html. AP. PWM, PDM, PS, PWetQ, 
PDQ, PCQ. PWarmQ 数据 来 自 1961—1990 4E, x 
[i] 4p XE XE Jy 30"(— 1 km), 来 源 于 http://www. 
worldclim.org. REH, WD 数据 来 自 1961—1990 年 ， 
空间 分 辨 率 为 10'， 来 源 于 http:/www.geodata.cn/ 
Portal. 

生境 异 质 性 假说 常用 生境 类 型 或 者 海拔 指标 
衡量 (Currie, 1991; Huang et al, 2011; Kerr & Packer, 
1997). KARHE PERS] T8848 352870. (land 
use and land cover change, LUCC)、 平均 海拔 (mean 
elevation, ELEmean)> 5 iE if d (maximum elevation, 
ELEmx). IRMA (minimum elevation, ELEmin) o 
拔 数据 来 源 http://eros.usgs.gov, 空间 分 辩 率 为 30'。 
植被 类 型 数据 来 源 国家 自然 科学 基金 委员 会 “中 国 
西部 环境 与 生态 科学 数据 中 心 ”http://westde. 
westgis.ac.cn， 空 间 分 辨 率 为 1 km?。 
133 统计 分 析 

不 同 空间 尺度 或 面积 对 物种 分 布 格局 的 形成 有 
重要 影响 (He & Legendre, 1996; Willis & Whittaker, 
2002)。 为 了 消除 空间 尺度 或 面积 对 物种 分 布 格局 的 
影响 , 本 研究 将 中 国 国土 划分 为 100 km 100 km 的 
栅 格 系统 (Huang et al, 2011); 为 了 防止 海岸 线 和 边 
界 地 区 不 完整 栅 格 影响 后 续 统计 分 析 ,， 去 除了 不 足 
完整 栅 格 面积 75% 的 栅 格 (Huang et al, 2011); 为 了 
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消除 岛屿 的 影响 (Connor & McCoy, 2001; Schall & 














N25?5ll N47°, 中央 经 度 E105°)。 同 时 , 所 有 生态 因 



































Pianka, 1978)， 去 除 所 有 岛屿 栅 格 。 本 研究 利用 
ArcGIS 软件 将 蛇 类 物种 分 布 数据 生成 每 个 物种 的 
Arc/Info 栅 格 图 有 某 一 物种 分 布 的 栅 格 赋值 为 “1?， 
无 某 一 物种 分 布 的 栅 格 赋值 为 “0”， 然后 登 加 所 有 物 
种 的 栅 格 图 以 生成 物种 丰富 度 格局 ， 最 后 投影 为 
Alberts Equal-Area Conic Projection ( 双 标 准 纬 线 : 



























































子 投影 重 采 样 为 100 kmX 100 km 并 提取 数据 。 

由 于 生态 因子 间 的 强 相 关 性 会 给 后 续 多 元 回归 
分 析 带 来 偏差 因此 研究 中 尽 可 能 选择 相关 系数 低 
的 变量 并 且 每 个 假说 至 少 有 一 个 变量 代表 (Huang et 
al 2011)。 于 是 , 在 SPSS 软件 中 对 生态 因子 进行 两 两 
之 间 相 关 分 析 后 , 选择 相关 系数 <0.7 的 变量 ( 表 2). 







































































表 2， 生 态 因子 的 相关 矩阵 (r <0.7) 


Tab.2 Correlation matrix of ecological factors (r <0.7) 





ELEmin WIND ISO TAR PS PCQ NDVI 

ELEmin 1 0.137** 0.611** -0.089** 0431** -0.390** -0.555** 

WIND 1 -0.224** 0.497** 0.346** -0.284** -0.483** 

ISO -0.518** 0.283** -0.188** -0.112** 

TAR 1 0.357** -0.546** -0.510** 

PS 1 -0.675** -0.406** 

PCQ 1 0.644** 
NDVI 1 





**; 相关 性 在 0.01 的 水 平 显 著 ( 双 尾 )。Correlation is significant at the 0.01 level (2-tailed). 
*#: 相关 性 在 0.05 的 水 平 显著 ( 双 尾 )。Correlation is significant at the 0.05 level (2-tailed). 











统计 分 析 在 SAM v4.0 (Spatial Analysis in 
Macroecology, http://www.ecoevol.ufg.br/sam) 
(Rangel et al, 2006, 2010) 软 件 中 完成 。 为 了 消除 回 
归 残 差 空间 自 相 关 ,， 本 研究 采用 主轴 邻 距 法 (PCNML 
principal coordinates neighbor matrices) 获 得 基于 特 
征 值 的 空域 数据 (eigenvector-based filters) 并 与 生态 
因子 一 起 进行 多 元 回归 分 析 ， 该 方法 能 有 效 降低 回 
归 残 差 空间 自 相 关 (Diniz-Filho & Bini, 2005)。 为 降 
低空 间 自 相关 性 ,本 文选 择 最 小 残 差 的 Moran's I 
指数 的 空域 文件 (filters)(Diniz-Filho & Bini, 2005), 
将 其 与 生态 因子 组 成 多 元 回归 模型 .Moran's 7 指数 
取 值 范围 为 -1 一 1, *-1" 意味 着 强烈 的 负 空间 自 相 
关 ，“0” 意 味 着 无 相关 性 , “1” 意 味 存在 强烈 的 































































































































































































正 空间 自 相 关 。 该 多 元 回归 分 析 包 括 了 空间 因素 , r 
由 三 部 分 组 成 : 仅 生 态 因子 (predictors only)、 仅 空 


























域 (filters only)、 两 者 (生态 因子 和 空域 )predictors+ 
filters)。 
为 了 探讨 5 个 假说 对 中 国 蛇 类 物种 丰富 度 分 布 




















的 模型 ， 是 在 模型 构建 中 偏差 和 精度 之 间 最 优化 的 
权衡 (Burnham & Anderson, 1998)。 当 所 有 模型 的 
AICc 值 都 很 接近 时 ,用 候选 模型 和 最 优 模型 之 差 
(AAICc) 衡 量 相关 模型 拟 合 度 。AAICc 越 小 ,候选 模 
型 越 接 近 最 优 模 型 。 即 ， 当 AAICc 科 2 时 ， 数 据 强力 
支持 该 模型 ， 当 AXAAICOXT 时 , 数据 一 定 程度 支 
持 该 模型 ; 当 AAICc>10 时 ,数据 不 支持 该 模型 。 
使 用 赤 池 权重 (Akaike weights, wi) 估 计 的 模型 平均 
数 来 解决 模型 选择 的 不 确定 性 (Bumham & 
Anderson, 1998). 


2 fd R 


2.4 物种 丰富 度 地 理 格局 

由 图 2 可 见 中 国 蛇 类 物种 丰富 度 热点 地 区 在 横 
断 山脉 南部 、 南 岭 、 武 夷 山脉 、 大 类 山 - 巫 山 山脉 、 
青藏 高 原 东 南 缘 、 秦 岭 - 大 巴山 脉 、 藏 南河 谷 ( 雅 
藏 布 江 下 游 地 区 )， 以 及 北方 的 长 白山 脉 、 小 兴安 岭 
和 西部 的 天 山 山脉 西部 等 。 中 国 蛇 类 物种 丰富 度 并 































































































































































































格局 的 解释 程度 ， 本 研究 对 各 假说 包含 的 所 有 变量 
采用 主轴 邻 距 法 分 别 进行 多 元 回归 分 析 ， 同 时 另 增 
加 两 个 混合 模型 : 大 混合 模型 (所 有 变量 )， 小 混合 
模型 (多 元 回归 分 析 中 P<0.05 的 变量 )。 采 用 样本 大 
小 矫正 过 的 赤 池 信息 准则 (AICc) 来 衡量 模型 拟 合 度 
优良 性 。 最 小 AICc 值 的 模型 是 最 简约 或 优先 考虑 






























































不 随 纬度 升 高 而 单调 降低 ,而 是 在 N22。 一 28" 为 最 
EE, 在 N28°~37° 随 纬度 升 高 而 显著 降低 ,在 ~ 
N37°* 有 一 个 低谷 (图 3A); 在 E75?—90? RAKK, 
在 E90°~~110° 随 着 经 度 升 高 而 上 升 ， 在 E110?— 
128° 随 着 经 度 升 高 而 降低 ， 在 E128? 之 后 随 着 经 度 
升 高 而 上 升 ， 丰 富 度 最 高 值 位 于 一 E110" 处 (图 
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图 2 ”中国 蛇 类 物种 丰富 度 栅 格 图 (100 km X 100 km) 


Fig.2 Grid map of snake species richness patterns in China at 100 km spatial resolution 


3B)。 结 合 图 2 可 知 中 国 蛇 类 物种 丰富 度 最 高 的 地 
区 出 现 东 洋 界 亚热带 、 热 带 地 区 , 物种 丰富 度 较 低 
的 地 区 在 青藏 高 原 、 北 方 草原 荒漠 地 区 、 黄 淮 平原 、 
两 湖 平原 、 好 阳 湖 平原 等 地 。 
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图 3 中 国 蛇 类 物种 丰富 度 与 纬度 (A)、 经 度 (B) 
关系 (100 km X 100 km) 

Fig.3 Relationship of snake richness (100 kmx100 km) in 
China with (A) latitude and (B) longitude 
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2.2 ”物种 丰富 度 地 理 格 局 与 生态 因子 的 关系 

主轴 邻 距 法 分 析 的 截 尾 距离 (truncation distance) 
为 142.265 km。 选 择 最 小 残 差 的 Moran's 了 指数 的 
衬 域 文件 将 其 与 生态 因子 组 成 多 元 回归 模型 ， 有 效 
消除 了 空间 自 相 关 性 (图 4)。 
























* 物种 丰富 度 Species richness 
* 回归 残 差 Residuals of regression 


500 1000 
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类 间距 离 Distance class 











VRI 


4 蛇 类 物种 丰富 度 和 多 元 回归 残 差 的 Moran's 7 指数 图 

Fig. 4 Moran’s J correlogram for snake species richness 
and the residuals of multiple regression with 
ecological predictors 


使 用 空间 特征 值 模型 的 多 元 回归 分 析 能 够 解 
释 蛇 类 物种 丰富 度 总 体 56.5% I] 28 45 (72—0.565; 
F-66.86; P=0), ESAT MERE T 39.9%, IXE ER 
解释 了 8.6%， 两 者 (生态 因子 和 空域 ) 共 同 解释 了 
8.1%。 在 所 有 生态 因子 中 ， 归 一 化 植被 指数 和 最 冷 
季 降 水 量 是 影响 中 国 蛇 类 物种 丰富 度 最 大 的 因子 ， 
且 它 们 均 与 中 国 蛇 类 物种 丰富 度 呈 正 相 关 , 年 温差 
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与 生态 因子 的 关系 
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是 影响 中 国 蛇 类 物种 丰富 度 次 大 的 因子 , 与 中 国 蛇 C 
3.4 ”中国 蛇 类 物种 丰富 度 地 理 格局 








类 物种 丰富 度 呈 负 相 关 ( 表 3)。 





2.3 生态 假说 对 物种 丰富 度 地 理 格 局 的 解释 
基于 亦 池 信 息 量 准则 的 模型 选择 方法 ， 小 混合 

















模型 (多 元 回归 分 析 中 P«0.05 的 变量 ) 的 AAICc 最 低 ， 

















赤 池 权重 为 1036 4), 是 解释 中 国 蛇 类 物种 丰富 度 格 

















局 的 最 优 模型 。 其 他 模型 (如 : 大 混合 模型 





、 生 境 异 


质 性 假说 、 生 产 力 假 说 、 环 境 能 量 假说 、 寒 冷 妨 而 
假说 、 水 能 动态 假说 ) AAICc 均 >10， 赤 池 权 重 也 


很 低 。 


表 3 ”中国 蛇 类 物种 丰富 度 与 生态 因子 的 多 元 回归 分 析 


Tab.3 Snake richness regressed against environmental 


variables in China 
生态 因子 











Ecological factors b : Bá 
最 冷 季 降水 量 PCQ 0.098 8.915 <0.001 
归 一 化 植被 指数 NDVI 10.934 3.345 <0.001 
年 温差 TAR -0.171 -2.557 0.011 
季节 性 降水 量 PS 0.027 1.437 0.151 
最 低 海拔 ELEmin <0.001 —0.907 0.365 
风速 WIND -0.161 —0.385 0.700 
等 温 性 ISO -0.015 —0.216 0.829 











b: 直线 回归 方程 系数 。The linear regression coefficients. 
未 显示 基于 主轴 邻 距 法 而 获得 空域 数据 值 (Filters obtained by Principal 
Coordinates Neighbor Matrices are not shown) (270.565, Fl(17.874)=66.86, 








P-0). 











t 直线 回归 t 检 验 值 。The linear regression t-test values. 
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2000)， 蛇 














岂 类 也 不 例外 。 中 国 





高 的 地 











种 丰富 度 具 有 空间 异 质 性 特征 (Gaston， 
区 域 蛇 类 物种 丰富 度 最 


热带 地 区 ， 如 : 武夷 山 


















































脉 、 横 断 山 脉 南 部 、 南 岭 、 藏 南河 谷 (雅鲁藏布江 下 





游 地 区 )、 大 关山 -巫山 山脉 、 





岭 - 大 巴山 脉 (图 2)。 这 些 地 




















RRRA 


Jie, 天 山 山脉 西部 等 地 ( 
同 特点 也 是 山区 , 与 北方 地 


E fg a AER 








青藏 高 原 东 南 缘 和 秦 

















区 均 为 南方 山区 , 与 中 
上 量 、 归 一 化 植被 指数 高 的 地 区 (图 5) 





。 男 外 北方 丰富 度 较 高 地 区 为 长 白山 脉 、 
图 2), 这 些 地 区 的 共 
区 归 一 化 植被 指数 、 引 








均 温 高 





Acus px 
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的 地 区 (图 DEAR. T 





藏 高 原 、 北方 草原 











、 大 兴安 岭 地 区 是 ! 


呢 类 物种 丰富 度 低 的 地 

















pe, P 
被 指数 、 


降水 量 、 归 一 化 植被 指数 与 蛇 类 物种 
关 。 最 冷 季 降 水 量 高 能 有 效 保 记 














省 严寒 或 干旱 ， 与 最 冷 季 降 雨量 








、 归 一 化 植 








年 均 温 低 的 地 区 ( 医 











5) EUR. 最 冷 季 























眼 期 间 有 足够 的 水 分 支撑 蛇 类 











保证 蛇 类 冬眠 所 需 环境 湿度 ， 帮 





























IRH; 归 一 化 植被 指数 高 能 
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:这 些 地 区 在 蛇 类 冬 











捕食 对 象 的 生存 ， 


X uk 








Ji yi NUR E SE AC 


上 促进 蛇 类 食物 链 下 端的 


动物 大 量 生长 繁殖 ， 增 加 蛇 类 数量 和 种 数 。 年 温差 








则 与 蛇 类 物种 丰富 度 呈 负 相 关 。 
效 的 体温 调节 机 制 , 但 作为 冷 


















































气温 范围 
生长 发 育 (Sun, 2001), 


得 蛇 类 








表 4 ”模型 选择 (K 表示 模型 变量 数 , 不 包括 截 距 ; wi 表示 权重 ) 
Tab.4 Summary of the model selection approach based on AICc w; was the Akaike weight (K was the number 
of variables of the model, not including intercept) 





























虽然 冬眠 是 蛇 类 有 
动物 ， 蛇 类 的 最 适 


是 20—30 '"C(Zhao, 2006), 过 大 温差 会 阻 
因此 年 温差 高 的 地 区 








生态 假说 











Ecological hypothesis h OE Bees K wi 
H 型 
大 混合 模型 0.57 6085.32 40.89 31 0.00 
Big mixed model 
HE 
。 EBOURTEBUR 0.56 6084.68 40.26 4 0.00 
Habitat heterogeneity hypothesis 
Hue 
生产 力 假说 0.56 6069.74 25.31 2 0.00 
Productivity hypothesis 
境 能 量 
环境 能 量 假说 0.56 6069.61 25.18 4 0.00 
Ambient energy hypothesis 
3e» Zi 
耐候 说 0.57 6065.40 20.97 10 0.00 
Freeing A hypothesis 
bbh 
水 能 动态 假说 ; 0.56 6064.55 20.12 13 0.00 
Water-energy dynamic hypothesis 
小 混合 模型 
Small mixed model based on the significant 0.57 6044.43 0.00 3 1.00 


variables (P«0.05) from the multiple regression 
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图 5 中 国 部 分 生态 因子 分 布 图 
Fig. 5 Maps of some ecological factors in China 

A) 海拔 (ELE)，B) 最 冷 季 均 温 PCQ)，C) 归 一 化 植被 指数 NDVD; D) 年 温差 (TAR)。 

A) Elevation (ELE); B) Precipitation of coldest quarter (PCQ); C) Normalized difference vegetation index (NDVJ); D) Temperature annual range (TAR). 





























蛇 类 物种 丰富 度 低 , 年 温差 低 的 地 区 蛇 类 物种 丰富 一 致 (Huang et al, 2011), 但 也 存在 细微 差异 , 较 明 
度 高 。 显 的 地 区 在 北方 草原 苑 漠 和 天 山 山脉 西部 地 区 。 中 
采用 主轴 邻 距 法 的 多 元 回归 分 析 结 果 与 海拔 、 国 蜥 蝎 类 物种 丰富 度 格 局 与 年 均 霜 期 呈 负 相关 ,与 
最 冷 季 降水 量 分 布 图 、 年 温差 分 布 图 具有 良好 的 对 ， 海拔、 植被 类 型 数 及 潮湿 天 数 频率 呈正 相关 (Huang 
应 关系 , 但 局 部 区 域 除外 ， 如: 东洋 界 的 黄 淮 平原 、 et al, 2011)， 而 年 均 霜 期 、 植 被 类 型 数 及 潮湿 天 数 频 
两 湖 平 原 、 郡 阳 湖 平原 等 区 域 昵 类 物种 丰富 度 很 率 与 本 研究 中 所 选 的 蛇 类 各 生态 因子 相关 性 很 高 ， 
低 。 这 些 区 域 的 晰 蝎 类 、 哺 乳 类 物种 丰富 度 也 低 (Lin ”海拔 相关 性 相对 较 低 ( 表 2), 海拔 变 幅 对 蜥 蝎 影 响 
et al, 2009; Huang et al, 2011)。 原因 可 能 是 这 些 区 域 。 显著 而 对 蛇 类 不 显著 , 与 Shen (1985) 的 结论 一 致 。 
历史 上 受 河 患 影响 严重 且 农 业 发 达 (Gong, 2002; 。 ”这 可 能 是 蛇 类 、 蜥 蝎 类 在 这 些 区 域 的 物种 丰富 度 存 
Lin & Mo, 2011; Zhao et al, 2007)。 中 国 扑 行动 物 濒 。 ”在 微小 差异 的 原因 之 一 ; 其 次 可 能 是 由 于 两 类 群 本 
危 的 主要 原因 是 食用 和 医用 造成 过 度 利用 ， 乱 采 滥 。”” 底 资料 调查 不 一 致 所 致 ， 因 而 有 待 进一步 研究 。 
伐 导 致 生境 破碎 化 , 污水 和 农药 污染 (Li & Wllcove， — 3.2. ”中国 蛇 类 物种 丰富 度 地 理 格局 的 成 因 及 其 复 































































































































































































2005)。 因 此 这 些 地 区 蛇 类 物种 丰富 度 低 , 除了 河 患 杂 性 
导致 外 , 还 可 能 是 由 人 类 影响 所 臻 。 具 体 原因 有 竺 生态 假说 中 , 除去 生境 异 质 性 假说 之 外 的 假说 
进一步 研究 。 均 属于 能 量 假说 (energy hypothesis), 该 假说 认为 气 


中 国 蛇 类 物种 丰富 度 分 布 格局 与 蜥 蝎 类 较为 候 的 能 量变 化 在 初期 产生 了 全 球 丰 富 度 梯 度 ， 能 量 











ps 


4 期 3k 波 等 : 上 















































2003a; von Humboldt, 1808). Hawkins et al (2003a) 
研究 得 出 能 量 组 成 的 气候 模型 能 更 好 地 解释 大 尺 
度 丰 富 度 格 局 , Guisan & Hofer (2003) 也 认为 中 等 尺 
度 能 量 组 成 的 气候 模型 比 地 形 模型 (地 形 为 生境 异 
质 性 假说 的 成 分 之 一 ) 更 能 解释 息 行 动物 物种 丰富 
度 格局 ,而 生境 异 质 性 假说 很 难 在 大 尺度 上 解释 和 候 
行动 物 空间 分 布 格局 (Kiester, 1971; Rodriguez et al, 
2005)。 但 也 存在 其 他 观点 。 中国 扑 行 动物 物种 丰富 
度 格局 成 因 研究 中 ,有 的 既 支 持 生境 异 质 性 假说 又 
支持 能 量 假说 (Huang et al, 2011; Qian et al 2007; 
Shen, 1985), A IJ p Sc fr Be Œ fix Vi (Fu. et al, 2007; 
Qian et al, 2007; Wang, 2004). 

中 国 蛇 类 物种 丰富 度 格 局 的 强 影 响 因子 为 最 
冷 季 降水 量 和 归 一 化 植被 指数 。 最 冷 季 降水 量 与 蛇 
类 物种 丰富 度 分 布 格局 呈正 相关 , 该 研究 结果 支持 
水 能 动态 假说 。 美 洲 的 珊瑚 蛇 属 (Terribile & 
Diniz-Filho, 2009) fI "F Ed f £T 2/) 4] ( Fu et al, 2007; 










































































































































































可 利用 性 导致 并 维持 物种 丰富 度 梯 度 (Hawkins et al, 
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量 。 这 说 明 物 种 丰富 度 格 局 影响 因子 多 样 ， 而 现 有 
J 假说 解释 度 低 。 如 果 进 行 赤 池 信息 量 准 则 的 模型 
选择 时 去 除 小 混合 模型 ， 得 出 的 解释 模型 将 依次 是 
水 能 动态 假说 、 赛 冷 忍 耐 假说 、 环 境 能 量 假 说 及 生 
产 力 假说 。 出 现 这 一 情况 的 原因 就 是 各 假说 中 生态 
因子 成 分 的 影响 。 很 多 生态 因子 之 间 关 系 密切 ,各 
假说 之 间 也 紧密 联系 (Allen et al, 2007; Wang et al, 
2009)， 各 假说 或 模型 受 生 态 因子 影响 很 大 ， 稍 有 不 
慎 则 结果 巡 异 ,如 本 研究 中 的 大 、 小 混合 模型 。 因 
此 , 应 当 重 视 各 假说 所 包含 影响 因子 的 选择 。 建 议 
先 参考 前 人 的 研究 , 选择 相关 假说 生态 因子 ， 尽 可 
能 全 面 获取 生态 因子 数据 及 后 各 个 生态 因子 的 变 
量 (如 最 低 值 、 最 高 值 、 平 均值 、 变 异 系数 等 )， 在 统 
计 分 析 前 根据 需要 去 掉 相 关 性 高 的 因子 (如 本 文 去 
挥 相关 性 >0.7 的 因子 )， 以 免 影响 统计 结果 。 

其 次 ,尺度 复杂 性 。 丰 富 度 格局 会 随 尺 度 变 化 ， 
局 域 尺 度 进程 也 影响 区 域 尺 度 丰 富 度 的 观测 
(Gaston, 2000)。 尺 度 具 有 变异 性 , 各 生态 因子 在 不 
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Qian et al, 2007; Wang, 2004) 相 关 研 究 也 支持 该 假 
说 。 该 假说 主要 基于 水 分 对 生物 体能 量 利 用 过 程 的 
影响 (Wang et al, 2009), 能 有 效 解 释 整 个 地 球 生物 
多 样 性 梯度 (Hawkins et al, 2003a), 但 对 动物 的 影响 
可 能 是 通过 对 植物 生产 力 和 生物 量 的 影响 而 间接 

































































同 尺度 表现 出 不 同 的 特征 (Fu et al, 2007; Gaston, 
2000; Lü & Fu, 2001)， 即 对 物种 多 样 性 分 布 格局 的 
文 持 度 ， 也 影响 各 假说 对 分 布 格局 的 解释 (Hawkins 
et al, 2003a). Hu et al (2005) 分 析 了 非 环境 类 因子 对 
物种 数目 协同 变化 的 影响 发 现 , 不 同 鸟 兽 间 的 协同 
































实现 (Hawkins et al, 2003a)。 归 一 化 植被 指数 对 蛇 类 
丰富 度 格局 的 影响 则 一 定 程 度 上 证 明了 这 种 间接 
性 。NDVI 是 生产 力 假说 的 重要 内 容 ， 本 研究 表明 
































变化 不 存在 面积 效应 , 但 其 空间 尺度 间 的 相关 性 差 
异 明显 。 本 研究 与 后 续 各 尺度 的 研究 (结果 男 文 发 表 ) 
表明 在 小 尺度 与 较 大 尺度 上 物种 丰富 度 格局 影响 






































NDVI 在 一 定 程度 支持 生产 力 假说 , 通过 控制 区 域 
内 的 蛇 类 食物 数量 来 控制 蛇 类 种 群 密度 ,进而 影响 
其 物种 丰富 度 水 平 。 而 研究 中 代表 生境 异 质 性 假说 
的 土地 履 盖 类 型 数 、 平 均 海 拔 、 最 高 海拔 、 最 低 海 
拔 与 蛇 类 物种 丰富 度 格局 相关 性 很 小 ， 因 此 本 研究 
文 持 能 量 假说 , 不 文 持 生境 异 质 性 假说 。 

物种 丰富 度 分 布 格局 成 因 研 究 中 虽然 出 现 了 
多 种 假说 或 者 理论 , 但 仍然 没有 单一 假说 能 充分 解 
释 (Gaston, 2000; Cox & Moore, 2010)。 本 研究 表明 









































































































































因子 不 同 但 存在 关联 ， 局 部 尺度 影响 因子 也 会 对 区 
域 尺度 产生 影响 (图 5),。 因 此 将 来 需要 研究 不 同 尺度 
对 物种 丰富 度 格局 的 影响 ， 以 及 不 同 尺 度 上 各 假说 
对 物种 丰富 度 格 局 成 因 的 解释 。 

最 后 ， 人 类 因素 影响 的 复杂 性 。 人 类 因素 影响 
属于 生物 因子 影响 之 一 ， 从 生物 干扰 、 食 物 竞争 、 
捕食 等 作用 上 直接 或 间接 影响 物种 分 布 格局 。Li & 
Wllcove (2005) 研 究 得 出 资源 过 度 利用 、 "EE CE 
以 及 各 种 污染 是 我 国 两 栖 扑 行动 物 濒危 的 主要 原 












































































































































各 种 假说 并 不 相互 排斥 ， 而 代表 各 假说 的 生态 因子 
共同 影响 中 国 蛇 类 物种 丰富 度 格局 。 本 研究 使 用 空 
间 特 征 值 模型 的 多 元 回归 分 析 能 够 解释 中 国 蛇 类 
物种 丰富 度 总 体 56.5% 的 变化 ， 尚 有 43.5% 不 能 解 
释 。 其 复杂 性 主要 表现 在 以 下 三 个 方面 。 

首先 , 各 假说 中 生态 因子 的 复杂 性 。 众 多 模型 
中 ,小 混合 模型 最 能 解释 中 国 蛇 类 物种 丰富 格局 ， 
因为 它 包括 了 在 多 元 回归 分 析 中 所 有 P<0.05 的 变 































































































| 









































因 。 另 外 ， 人 类 生产 生活 导致 的 气候 变化 将 影响 动 
物 的 行为 和 物候 (Bull & Vogt, 1979; Gargaglioni & 
Milsom, 2007)， 改 变 分 布 范 围 ， 继 而 改变 丰富 度 其 
至 导致 物种 灭绝 (Pounds et al, 1999; Sun & Zhang, 
2000)。 蛇 类 具有 重要 的 经 济 、 科 研 价 值 ， 尽管 在 
2000 年 出 版 《国家 保护 的 有 益 的 或 者 有 重要 经 济 、 
科学 研究 价值 的 陆 生 野生 动物 名 录 》 所 列 的 中 国 
395 PRETA, 蛇 类 就 有 207 种 ， 是 当时 我 国有 
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分 布 记录 的 所 有 蛇 类 , 但 我 国 





蛇 类 受 威胁 程度 很 高 ， 














几乎 所 有 蛇 的 整体 和 器 官 都 能 作为 中 药 、 工 艺 、 食 
品 材 料 (Zhao, 2006)。 因 此 ,必须 重视 人 类 因素 对 





我 国 蛇 类 丰富 度 格局 的 影响 。 

中 国 蛇 类 物种 丰富 度 地 理 格局 成 因 研 究 有 
其 复杂 性 ， 因 此 ， 在 进一步 研究 中 需要 重视 各 假 
说 中 影响 因子 的 选择 和 人 类 活动 的 影响 ,并 在 不 同 
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空间 尺度 上 对 蛇 类 物种 丰富 度 地 理 格局 进行 综合 
分 析 。 
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